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a 
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ABSTRACT 


An  experimental  investigation  has  been  made  to  determine  the 
pressure  distribution,  shock  shape,  and  shock  detachment  distance 
for  spheres  and  the  latter  two  characteristics  for  flat»faced  bodies 
in  a  heated  argon  flow  where  =  4  to  14,  Tq  =  1900  to  4100“K,  and 
Re2  =  ^5  to  225.  *> 

The  modified  Newtonian  approximation  for  the  pressure  distribu¬ 
tion,  which  is  strictly  an  empirical  relationship,  gives  good  results 
when  applied  to  the  first  60 -deg  of  a  hemisphere  under  conditFons  where 
the  body  boundary  la,yetf  an«f  the  shogk  layer  merge. 

The  natura?  fte'^tf  visdallzaUolS  produced  as  a  result  of  the  high  total 
temperatur®  an^  COnseqtieftt  eycitalien  of  the  argon  enabled  this  study  to 
be  Wade,  Photographs  fahen  of  She  shocks  generated  by  a  series  of 
spheres  and  tlat*|aced  bodies  Wer®  analysed  with  a  pfioto- densitometer 
to  determine  the  shofib  ehape  and  shook  detafthinent  distance®  The  blast 
analogy  predicts  a  difference  In  the  shape  e>f  shoOks  in  argon  and  air  at 
high  Mach  numbers  ajjd  numbers.  An  emfprical  relationship 

proposed  by  bove  was  used  to  calculate  the  shock  shape  iij  air  for  Mach 
numbers  correspojiding  to  those  t>£  the  present  Investigation. 

When  this  shape  was  compared  to  tne  shape  measured  in  ar^n  at 
the  same  Mach  number,  it  was  found  that  the  difference  in  sjlape  was  in 
agreement  with  that  predicted  by  the  blast  analogy,  i.  e. ,  iog  a  {onstant 
value  of  x/D  then  ty/OJargoii/<//C))air  *  1.  1^.  Also,  it  appears  even  for 
the  lowest  Reynolds  number  tested  that  nondimensionalized  shock  shape 
is  independent  of  Reynolds  numher  for  spheres  and  flat-faced  bodies  when 
the  effect  Of  Reynolds  numbei*  on  detachment  distance  is  taken  into  account. 

A  study  of  the  shock  detaclunenl  distance  indicates  that  the  boundary 
and  shock  layers  were  either  inclplently  or  fully  merged  In  these  experi¬ 
ments.  FurthercBore,  for  the  Reynolds  number  range  of  the  present 
tests,  the  detachment  distance  is  a  ftiaetion  of  wall  temperature,  Reynolds 
number,  Mach  number,  and  body  shape.  As  the  Reynolds  number  de¬ 
creases,  the  shock  detacbmesit  distance  increased  to  valPes  more  than 
double  the  "Inviscid"  values.  Also,  for  the  Water-cooled  bodies,  there  is 
evidence  to  confirm  the  decrease  iij  detachment  distance  to  less  than  the 
inviscicj  value,  as  predicted  by  Ho  aijd  Probslejn,  before  the  increase  as 
mentioned  above  for  the  lowest  Reynolds  numbers. 
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NOMENCLATURE 


o 


A,  a,  and  n 

Constants  in  Eq.  (6) 

Nose,  draf^  coefficient 

% 

TJ  fi*"  i  ^  ~  ^Od) 

Pressure  coeflicient  =  - — 

D  U  ^ 

D 

Body  diarPLoter 

s 

Sho<<®?:owaye  thickness 

Jo 

o 

^$ast  ana’ogy  Skinstar.t.  Eq.  (7) 

0  0 

% 

PaotoiOg  Eq.  ^1) 

Factor,  Eq.  (5? 

M2 

Magh  numlger  ?n;nnediate9y  behind  a  n9rnial 

Free-stream  Mach  number 

1% 

Mass  f ;  o  w  t  e 

o 

Pressure 

P2 

PressQ?re  ?mmed?ate?y  downctfeam  of  a  n^r^Sifcat  ^.ock 

P® 

Freeostream  pressure 

^li^ssure 

^°c 

b^eservoir  press<>re  -  ©old  flow 

R^ervofj*|)ressur^  -  hot  flow 

BOc^  ra(7?ifs  nos^  radios  of  curvalftre  for  sphere) 

Reyndjds  n#niber  bas^  on  body  diameter  and  5®n(?it?on 
iijpmediat^y  downstf^m  of  a  normal  shock®  For 
pujrp^es  of  (*a?i?ulating  Re2  this  shoij^  is  assumed  to  be 
ifnaffegted  by  iTie  boundary' ?ayer. 

5 

^olegular  spe^  ratio 

9 

Ve 

O 

Reservoir  temperature  ^  gas 

^esQi’vaffr  Temperature  of  gas  -  (gpld  fl^w^ 

Reservoir  temperature  of  gas  -  hot  flow 

TV 

Boc^  wal?  temperature 

Veloci^'  immediately  downstream  of  a  normal  shoigk 

0 

Free- stream  velod^itv 

€/  O 
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ix 


0/ 


AE»C.TDR.63-21 


X,  y 
Xb--  y 
y 

A 


s 

6*  6 
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A- 

o 

/^w 

P2 

Pm 

Pw 


o 


o 


Velocity  gradient  at  the  stagnation  point 

Ordinates  with  origin  at  the  apex  of  the  shock 

Ordinates  with  origin  at  the  apex  of  the  body 

Specific  heat  ratio 

Free- stream  specific  heat  ratio 

Shock  detachment  distance,  distance  from  the 
body  to  the  shock  leading  edge  measured  along 
the  axis  of  symmetry 

Boundary-layer  thi<®kness 

;@Qundary- layer  displa^ment  thickness 

Angle  between  the  local  normal  to  the  surface 
and  the  direction  of  the  ^mdisturbed  free-stream 
velocity 

Freeestrgam.  mean  free  path 

goeffTijfent  of  viscosity  basad^n  wa?^  temporaj^ra 
Density  ?rrtm^d?ately  downstrean^  of  a  norma?  sho^ 
me  -strean^densfl^ 

e  e 

Density  at  the  wa>? 
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1.0  INTRODUCTION 


The  shock  detachment  distance  in  front  of  blunt  bodies  at  high  Reyn¬ 
olds  numbers  has  been  studied  both  theoretically  and  experimentally. 

Van  Dyke  (Ref.  1)  presents  a  theoretical  analysis  for  a  series  of  perfect 
gases  having  specific  heat  ratios  of  1,  5/3,  and  7/5  and  a  Mach  number 
range  from  1«  2  to  infinity.  This  analysis  for  a  sphere  indicates  that  the 
shock  detachment  distance  is  primarily  a  function  of  the  density  ratio 
across  the  norma?  shock,  the  effect  of  specific  heat  variation  at  fixed 
density  ratio  ?s  small,  and  there  is  no  variation  with  Reynolds  number  in 
essentially  inviscid  flows.  In  this  analysis  it  is  implied  that  the  shock 
and  boundary- la^e^  thicknesses  are  small  compared  w?th  the  shock  de¬ 
tachment  d^stancQa  At  lov.'  Reynolds  numbers  the  shcfjK  and  boundary- 
layer  thicknesses  are  no  longer  small  ®pmpared  with  shock  detachment 
distance,  and  at  sufficiently  low  Reynolds  numbers  the  shock  and  bound¬ 
ary  layer  will  merge.  This  implfes  tha^ there  must  be  a  Reynolds  number 
below  whioh  the  shock  detachment  distant  becomes  a  function  of  Reyn¬ 
olds  numbe^as  v/e?!  as  density  ratio. 

Probstoin  ?.nd  Tf.---.-.-  V.  and  Levinsk^  a 

(Refs©  2-4)  have  made  theoretical  analyses  of  this  problem  in  the  incip¬ 
ient  merged *fayer  flow  regime.  These  analyses  have  one  thing  in  common: 
as  th^ Reynolds  number  is  reduced  through  this  r^ime  and  the  fully 
merged  regime  is  approached,  the  shock  deta^ment  distance  ultimately 
increases.  Reference  2  is  based  on  a  constant  density  assumption  for  the 
shock  layer.  Ref.  S  this  to  consideration  of  variable  density,  and 

Ref.  4  assumes  constant  density  but  differs  frorr^Re%  2  in  that  a  monatomic 
gas  (argon)  is  cons?d^ed. 

The  Low- Density,  Hypervelocily  (loDH)  Wi]:jd  Tunnel,  when  operated 
with  argon  as  the  worrt9hg  gas,  is  w^i  suited  foi:»?ow  Reynolds  number, 
shock  detachment  distanj^  investigations.  Wh^  th?fe  tiJnnel  is  operated 
•using  argonp  t?iere  is  a  natural  fi^w  visualization  thought  to  be  caused  by 
radiation  from  reiaxihg  metastable  argon  atoms.  Th?s  natural  glow 
enables  the  shock  shapes  in  front  of  bodies  to  be  seen  and  photographed. 
Moulic  (Ref,  5)  has  photographed  some  shock  shapes  in  a  low- density, 
supersonic  wind  tunnel  at  Berketey©  The  flow  visualization  in  that  case 
was  artifically  produced  by  infroduiyng  some  nitric  oxide  and  atomic 
oxygen  into  the  air  stream. 


Manusi^ipt  released  by  authors  January  19G3. 
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Also  of  fnterest  is  the  form  of  the  shock  wave  in  front  of  bodies  in 
low-density  flows.  There  is  a  considerable  amount  of  data  for  air  at 
high  Reynolds  numbers  but  very  little  for  other  gases  and  very  little 
in  any  case  for  low  Reynolds  numbers.  The  form  of  the  shock  wave  is 
predicted  reasonably  'well  by  the  blast  analogy  (Ref.  6),  and  it  was  be¬ 
lieved  to  be  of  interest  to  determine  whether  that  method  is  still  valid 
in  the  low-density  regime. 

o 

0 

For  high  Reynolds  numbers  the  pressure  distribution  over  a  sphere 
can  be  calculated  reasonably  accurately  by  using  the  modified  Newtonian 
approximation  (Ref.  7).  It  was  believed  interesting  to  ascertain  whether 
this  approximation  is  valid  in  a  flow  regime  where  the  shock  and 
boundary  layers  are  merging.  A  knowledge  of  the  form  of  the  pressure 
distribution  is  useful  because  the  stagnation  point  heat  transfer  and  the 
growth  of  ihe  laminar  boundary  layer  are  functions  of  the  velocity  gra¬ 
dient  at  the  stagnation  region  which  is  in  tui’n  a  function  of  the  pressure 
distribution. 


2.0  APPARATUS 

2.1  WIND  TUNNEL  DESCRIPTION  AND  PERFORMANCE 

The  LDH  Wind  Tunnel  is  a  continuous- type,  high-enthalpy  wind  tun¬ 
nel  in  operation  at  the  von  Karman  Gas  Dynamics  Facility  (VKF)  of  the 
Arnold  Engineering  Development  Center  (AEDC),  Air  Force  Systems 
Command  (AFSC),  U.  S.  Air  Force,  Briefly,  the  tunnel  consists  of  a 
d-c  arc  heater,  a  stftl  ing  chamber,  an  aerodynamic  nozzle,  a  test  cham¬ 
ber  with  instrumentatfon,  a  diffuser,  and  a  pumping  system.  A  sche¬ 
matic  drawing  identft';>'%g  som.e  of  the  components  is  show'n  in  Fig.  1 .  A 
complete  description  of  the  loDH  tunnel  is  given  in  Ref.  8. 

O  O 

In  the  pr  esent  tests  a  number  of  simiple,  water-cooled,  brass,  conical 
nozzles  with  a  30- deg  tctal  angle  were  used  to  accelerate  the  flow  to  super- 
son?!#  speeds.  The  brass  sections  have  throat  diameters  of  0.  200,  0.  397, 
and  0.  750  in.  whi^  produce  Mach  numbers  at  the  exit  plane  in  the  range 
from  4  to  8.  5  whei^ argon  is  the  med?bm.  They  can  be  used  in  conjunction 
w?th  an  aluminum  ^bne*(^us^im  extension  of  30  deg  total  angle  w'hich  con¬ 
tinues  the  expansfon  ®f  the  three  basic  nozzles  to  the  Mach  number  range 
®  from  7  to  14. 

o 

With  argon  as  the  working  gas  the  arc?  heater  operates  at  less  than 
15  kw.  Thc^gas  flow  rate  ranges  fromi  1  to  30  lb/ hr,  Tq  from  1900  to 
4100°K,  and  the  reservoir  pressure  ranges  from  1  to  6  psia.  At  the  nozzle 

®  o 

6 


2 


o 


o 


o 
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exit  and  some  distance  downstream,  there  is  a  clearly  visible,  light 
blue  jet  surrounded  in  most  cases  by  a  pink  region.  On  investigation 
with  an  impact-pressure  probe,  the  light  blue  region  was  found  to 
correspond  approximately  to  the  high-speed  core  of  flow.  The  colors 
in  the  flow  are  thought  to  be  caused  by  the  excitation  ot  argon  atoms  in 
a  metastable  state  since  the  enthalpy  is  too  low  for  significant  ionization. 
Good  in-flow  characteristics  to  the  aerodynamic  nozzles  under  test  were 
obtained  by  using  a  3-in.  -diam  stilling  chamber  upstream  of  the  nozzle 
throat,  having  a  length  of  4  in.  for  the  0.  200-  and  0.  397- in.  -diam 
throats  and  8  in.  for  the  0.  750-in.  -diam  throat. 

A  complete  description  of  the  gas  flow  control,  pressure  measuring 
system,  reservoir  temperature  estimate,  and  rotary  probe  holder  is 
given  in  Re%  S.  With  regard  to  the  reservoir  temperature  estimate,  the 
fact  that  argon  is  monatomic  considerably  simplifies  this  estimate  since 
perfect  gas  laws^pply.  Hence,  at  constant  mass  flow, 


2.2  TUHNEL  OPERATING  CONDITIONS 

The  usable  flow  regions  existing  with  the  aerodynamic  nozzles  were 
determined  by  making  impact-pressure  surveys.  These  nozzles  were 
operated  in  over-  and  underexpanded  slates,  and  for  all  the  flow  condi¬ 
tions  there  was  an  axial  Mach  number  gradient^  It  is  not  considered  that 
the  existence  of  an  axial  Mach  number  gradient  would  measurably  affect 
the  flow  conditions  at  the  front  of  a  model  provided  that  the  front  faces  of 
the  various  models  are  brought  to  the  same  axial  station  in  the  flow. 
Berndt  (Ref.  10)  has  made  some  est^iates  as  to  the  effect  of  conical  flow 
on  the  measured  shock  detachment  distance  for  high  Reynolds  number 
flows.  If  it  is  assumed  that  there  %  a  similar  effect  in  low  Reynolds  num¬ 
ber  flows,  then  Refjjp  10  ^ov/s  that  the  picsent  expei'imental  results  for 
shock  detachment  distance  cai^be  *9h  error  by  as  much  as  5  percent  on 
this  account. 


¥  PROCEDURE 


3.1  OPERATING  CONDITIONS 

The  desired  n*zz^e  flow  w^as  established  1^  adjusting  the  gas  flow- 
rate  and  arc-heaterflpower  input  to  bring  the  reservoir  conditions  to  the 
operating  level.  When  these  conditions  were  achieved,  tlie  model  under 
test  was  brought  to  a  predetermined  axial  station  on  the  nozzle  centerline. 
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3.2  PHOTOGRAPHIC  TECHNIQUE 

The  camera  used  in  this  series  of  tests  was  a  standard  4x5  Speed 
Graphic  with  fittings  such  that  a  Polaroid  back  could  be  attached  to  it. 

The  camera  was  mounted  outside  the  tunnel,  and  the  model  under  test 
was  viewed  through  one  of  the  windows  in  the  test  chamber.  A  13-in, 
focal  lengtii  lens  which  could  be  stopped  down  to  f/5.  6  was  used  with  the 
camera. 

The  relatively  low  intensity  of  the  glow  presented  some  problems 
in  defining  the  film  and  camera  settings  that  would  give  the  best  pictures. 

It  was  determined  by  means  of  a  trial  and  erior  procedure  that  good  re¬ 
sults  could  be  obtained  with  Kodalith  ASA  2  film.  It  was  found  that  with 
the  lens  stopped  down  to  f/5.  6  and  by  varying  the  exposure  from  two  to 
four  seconds,  dependent  on  the  model  position  with  respect  to  the  nozzle 
exit  plane,  good  quality  negatives  could  be  obtained.  A  typical  photo¬ 
graph  is  shown  in  Fi 

3.3  SPHERE  PRESSURE  MEASUREMENTS 

Because  of  the  small  size  of  the  model  (u.  5U-in,  diaml  and  the  need 
to  have  as  large  an  orifice  as  possible  from  pressure  lag-time  considera¬ 
tions,  it  was  not  practical  to  put  more  than  one  orifice  on  a  model  of  this 
size.  Therefore,  the  pressure  distribution  over  the  front  face  of  the 
sphere  was  obtained  by  using  a  number  of  models  each  with  one  orifice  at 
a  different  angle  to  the  horizontal  axis  of  the  sphere.  In  all,  eight  water- 
cooled  models  were  made  with  0.  052-in.  -diam  orifices  located  at  0,  15, 

30,  45,  60,  75,  90,  and  100  deg  to  the  horizontal  axis.  Each  of  the  models 
was  mounted  on  a  sting  which  could  be  attached  to  the  tunnel  probe  carrier. 
With  this  arrangement  it  was  possible  to  mount  the  models  such  that  the 
pressure  sensing  orifices  were  positioned  to  within  ±lo0  deg  of  the  nominal 
setting. 

For  each  flow  condition  the  front  face  of  the  sphere  was  brought  to  the 
predetermined  axial  position  and  the  pressure  was  measured  using  a  Pace 
transducer  with  an  eperating  range  of  ±0.  1  psi. 

O  -<.0  DISCUSSION  OF  RESULTS 

4.1  SHOCK  DETACHMENT  DISTANCE 

O 

4.1.1  Analysis  of  Photogrophs 

O 

Because  of  the  low-density  flow  conditions  existing  in  the  present 
tests,  it  is  no  longer  valid  to  assume  that  the  shock  wave  is  a  discrete 
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discontinuity.  For  example,  Levinsky  and  Yoshihara's  (Ref.  4)  anal¬ 
ysis,  in  the  region  of  the  stagnation  point,  indicates  that  at  the  lower 
Reynolds  numbers  the  variation  in  density  from  the  front  of  the  shock 
to  the  body  is  gradual  and  continuous.  If  it  is  assumed  that  the  density 
of  the  lilm  negative  in  this  region  is  proportional  to  the  gas  density, 
then  the  front  of  the  shock  will  be  determined  by  the  point  at  which  the 
film  density  starts  to  increase. 

A  photo- densitometer  was  used  to  determine  the  position  of  the 
shock  leading  edge,  the  front  face  of  the  sphere,  and,  where  possible, 
the  back  face  of  the  sphere.  The  shock  leading  edge  was  determined 
for  a  whole  series  of  off-axis  stations.  The  advantage  of  this  is  that 
it  provided  the  shock  shape  data  required  in  another  part  of  this  work, 
and  also,  by  draw'ing  a  smooth  curve  through  the  measured  points, 
errors  introduced  by  random  measuring  errors  in  stand-off  distance 
were  minimized.  This  type  of  error  is  particularly  critical  in  the 
determination  of  the  shock  detachment  distance  with  the  smallest  spheres. 
The  position  of  the  photo-densitometer  carriage  can  be  determined  to 
0,  f)01  cm;  however,  the  random  measuring  error  can  be  as  much  as 
+0.  005  cm,  which,  for  the  smallest  spheres,  implies  a  maximum  error 
of  approximately  ±10  percent  in  the  shock  detachment  distance.  This 
is  the  greater  error,  and  in  general  the  error  should  be  less  than  this 
for  a  negative  of  good  quality. 

4.1 .2  Sphere  Results 

Two  sets  of  spheres  were  studied  in  tliis  experiment:  one  set  con¬ 
sisted  of  uncooled,  stainless  steel,  ball  bearings  with  diameters  ranging 
from  0,  250  to  0.  G8G  in.  and  the  other  set  consisted  of  water-cooted 
brass  spheres  witli  diameters  ranging  from  0.  240  to  0.  749  in. 

A  preliminary  study  of  the  effect  of  wall  temperature  on  the  shock 
detachment  distance  indicated  that  the  detachment  distance  was  quite 
se^itive  to  relatively  small  clianges  in  wall  temperature.  Because  of 
this  temperature  sensitivity,  the  unc(>olcd  m.odels  were  allow'cd  to  reach 
their  equilibrium  temperature  before  a  photograph  was  taken.  The  time 
to  reach  thi^  equilibrium  temperature  was  determined  by  measuring  the 
sphere  temperature  with  a  chromei-alumel  thermocouple  attached  to  the 
rear  fa.ee  of  tne  sphere,  at  a  series  of  lime  increments  after  the  sphere 
was  brought  into  the  hot  stream.  It  was  found  from  this  investigation  that 
the  equilibrium  condition  was  attained  in  a  relatively  short  time  (  3  min.  ) 

and  that  the  wall  temperature  '.vas  approximately  1000°K.  Previous  ex¬ 
perience  with  j^faler-cooled  probes  in  the  LDH  Wind  Tunnel  has  shown  that 
the  water-cooled  spheres  should  quickly  reach  an  equilibrium  state  and 
tha'  the  wall  temperature  for  these  cases  is  approximately  300°K. 
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In  Figs.  3  and  4  the  results  of  the  determination  of  the  shock  de¬ 
tachment  distance  for  the  cooled  and  uncooled  spheres  are  plotted  in  the 
form  A/R)3  vs  Re2.  A  study  of  these  curves  shows  the  dependence  of 
shock  detachment  distance  on  Mach  number.  Reynolds  number,  and  wall 
temperature.  For  the  range  of  Reynolds  numbers  of  the  present  tests, 
the  shock  detachment  distance  increases  with  decrease  of  Reynolds  num¬ 
ber.  This  is  in  qualitative  agreement  with  the  theoretical  predictions  of 
Refs.  2-4  for  the  variation  of  shock  detachment  distance  with  Reynolds 
number  in  the  incipient  merged  layer  regime.  Furthermore,  it  is 
shown  that,  for  a  particular  Mach  number  and  Reynolds  number,  as  the 
v/all  temperature  increases  the  shock  detachment  distance  increases. 
This  is  in  qualitative  agreement  with  results  of  Ref.  4. 

Ii  is  of  interest  to  consider  the  relative  magnitu'^es  of  the  shock- 
layer  and  laminar  boundary- layer  thicknesses  at  the  flow  condition  of 
the  present  tests.  Ziering,  Ek,  and  Koch  (Ref.  11),  using  the  Boltz¬ 
mann  equation,  present  a  solution  to  the  problem  of  determining  the 
shock-wave  thickness.  A  comparison  of  experimental  data  obtained  on 
shock  thickness  in  argon  with  this  theory  indicates  that  it  gives  a  good 
estimate  of  the  shock-wave  thickness.  Reference  11  shows  that  the 
shock-wave  thickness  is  a  function  of  the  free-stream  Mach  number  and 
mean  free  path.  The  form  of  this  relationship  for  argon  is  shown  in 
Fig.  5.  Since  there  is  no  reason  to  doubt  the  validity  of  this  relation¬ 
ship  at  low  Reynolds  numbers,  the  shock-w'ave  thickness  has  been  de¬ 
rived  from  Fig.  5  for  the  range  of  flow  conditions  of  the  present  tests. 
Some  typical  values  are  shown  plotted  in  terms  of  body  radius  in 
Figs-  3  and  4.  It  can  be  seen  that  the  calculated  sliock-wave  thickness 
is  always  an  appreciable  percentage  of  the  measured  shock  detachment 
distance  at  these  low'  Reynolds  numbers,  ajid  hi  cases  it  is  siiown 

to  be  larger  ^han  the  measured  shock  detachment  distance.  This  may  be 
attributable  to  the  close  proximity  of  the  cooled  body  and  consequent  heat 
transfer  through  the  viscous  shock  layer.  This  possibility  is  not  con¬ 
sidered  in  the  theory  which  considers  isolated  shocks. 

The  laminar  boundary- layer  thickness  at  the  stagnation  point  can  be 
calculated  for  the  "thin"  boundary-layer  case  using  data  contained  in 
Ref.  12,  where  it  is  shown  that  the  laminar  boundary-layer  displacement 
thickness  at  the  stagnation  point  car.  Rewritten  in  the  following  form; 


k 


(1) 


where  K  is  a  factor  solely  dependent  on  the  ratio  of  wall  to  stagnation 
enthalpy  and  is  given  in  Ref.  12.  If  it  is  assumed  that  the  pressure  dis¬ 
tribution  in  the  region  of  the  stagnation  point  is  adequately  defined  by  the 
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modified  Newtonian  approximation,  then  the  velocity  gradient  at  the 
stagnation  point  for  a  hemisphere  can  be  written  in  the  form; 


,  r~. . 

1.7  -i-n~ 


‘  i>  ^2 

when  >  6.  0,  then  P2  >>  p,„  and  Eq.  (2)  may  be  rewritten  as 


U 


(2) 


(3) 


SubsliLuling  this  value  for  Ug  in  Eq.  (1),  and  writing  and  in  terms 
of  P2<  ^2-  ’i’2.  and  T^.  Eq.  (1)  becomes,  for  the  hemisphere. 


8" 


Rb 


(4) 


For  a  constant  value  of  y  the  product  of  the  teinio  containing  Tq/'Y^, 

M2j  and  P2/Po^  is  very  nearly  constant  for  a  wide  range  of  free- stream 

Mach  ??umbers  (5  <  <  20).  The  factor  K  is  a  function  of  the  ratio  of 

the  wall  to  stagnation  temperature,  and  K  changes  from  negative  to 

positive  as  T^/Tq  increases  from  zero  to  unity.  This  means  that  for 

highly  cooled  bodies  the  displacement  thickness  in  the  stagnation  region 

is  negative.  The  magnitude  of  the  displacement  thickness  is  a  function 

•  ,  V  I  S 

of  the  parameter  (  1', j  ^ 


In  the  present  context  it  is  considered  that  the  boundary-layer  thick¬ 
ness,  6,  is  a  more  meaningful  parameter  than  6*.  Referenc^  12  shows 
that  there  is  a  relationship  between  6  and  6*!'  which  is  a  function  of  wall  to 
stagnation  temperature  ratio  only.  Substituting  this  relationship  into 
Eq.  (4)  gives 


3  dli, 
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where  in  this  case  k  always  has  a  positive  value.  Typical  values  of  the 
boundary-layer  thickness  have  been  calculated  from  Eq.  (5)  and  are 
plotted  in  Figs.  3  and  4. 


For  the  same  Mach  number  and  Reynolds  number,  Eq.  (5)  predicts 
that  as  the  wall  to  stagnation  temperature  ratio  increases,  the  boundary- 
layer  thickness  decreases.  Equation  (4)  predicts  that  for  the  same  con¬ 
ditions  as  above,  the  displacement  thickness  increases  with  an  increase 
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in  wall  to  stagnation  temperature  ratio.  In  considering  the  case  where 
the  shock  and  boundary  layers  are  merged,  it  is  apparent  that  the  dis¬ 
placement  thickness  does  not  lose  its  significance;  but  the  total  thick¬ 
ness  probably  is  more  important  because  it  indicates  the  degree  of 
viscous  effect  in  the  shock  layer,  i.  e. ,  the  nearness  to  merging.  Nat¬ 
urally,  both  the  calculated  shock-wave  and  boundary-layer  thicknesses 
presented  here  are  not  exact  under  the  very  low  Reynolds  number  condi¬ 
tions  being  studied.  Tliey  are  presented  only  for  qualitative  comparisons. 

If  the  boundary-layer  thickness,  as  computed  by  a  theory  applicable 
to  the  case  6  «  R^,  is  considered  to  be  a  reasonable  approximation, 
then  Fig.  3  shows  that  all  of  the  results  of  this  test  are  in  the  fully 
merged  regime.  It  is  of  interest  to  note  that  the  form  of  the  measured 
shock  detachment  distance  variation  with  Reynolds  number  is  qualitatively 
similar  to  the  variation  of  shock-wave  and  boundary-layer  thickness  with 
Reynolds  number. 

In  the  merged-layer  cooled  model  case,  the  temperature  immedi¬ 
ately  behind  the  shock  is  reduced  to  a  value  less  than  the  stagnation 
temperature.  In  Ref.  11,  for  example,  shock- wave  thickness  is  com¬ 
puted  and  presented  in  terms  of  a  ratio  of  shock  thickness  to  free-stream 
mean  free  path.  Instead  of  referencing  shock-wave  thickness  to  the  free- 
stream  mean  free  path,  one  may  consider  a  value  of  mean  free  path  which 
is  the  average  of  the  free  stream  and  that  just  downstream  of  the  shock 
wave.  If  the  temperature  just  downstream  of  the  shock  is  reduced  be- 
cau.«e  of  wall  cooling,  then  the  average  value  of  the  mean  free  path  will 
be  reduced,  and  hence  the  shock  thickness  will  also  be  reduced.  This, 
plus  the  influence  of  cooling  on  density  in  the  shock  layer,  woiald  have  the 
effect  of  reducing  the  rate  of  increase  of  shock-layer  thickness  with  de¬ 
crease  of  Reynolds  number.  Some  confirmation  of  this  is  contained  in 
Figs.  3  and  4  when  the  data  for  M  6  are  compared.  There^fit  will  be 
noted  that  shock-layer  thickness  at  higher  Reynolds  numbers  is  greater 
for  the  uncooicd  body  which  theoretically  has  the  thinner  boundary  layer 
and  less  cooling  of  the  shock  layer. 

The  fact  that,  for  the  highest  Reynolds  numbers  tested  here,  the 
uncooled  sphere  shock  detachment  distance  is  greater  than  that  for  the 
cooled  spheres  may  be  explained  by  the  fact  that  in  the  first  case  the 
displacement  thickness  is  small  but  positive;  whereas  in  the  latter  case, 
it  is  small  but  negative--the  sum  of  these  two  values  accounting  for  the 
difference  in  measured  detachment  distance. 

A  method  of  correlating  the  present  results  with  those  obtained  for 
other  flow  conditions  is  suggested  in  Ref,  13  where  (A/Rb)/(A/R]3)j[jjvigcj[(j 
is  plotted  against  Re2-  Schwartz  and  Eckerman  (Ref.  14)  and  Vas,  et  al. 
(Ref.  15),  made  some  tests  to  determine  the  shock  detachment  distance 
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in  argon  and  helium,  respectively,  at  high  Reynolds  numbers,  the 
results  of  which  are  shown  in  b’ig.  6,  The  inviscid  values  for  (A/Rb) 
for  the  present  tests  were  taken  from  this  curve. 

Using  the  above-mentioned  parameter,  the  present  results  and  the 
theoretical  analyses  of  Refs.  2-4  are  compared  in  Fig.  7.  The  quali¬ 
tative  effect  of  wall  cooling  on  the  shock  detachment  distance  as  pre¬ 
dicted  by  Levinsky  and  Yoshihara  is  confirmed.  The  decrease  of  shock 
detachment  distance  to  a  value  less  than  the  "inviscid''  value,  as  pre- 
9  dieted  by  Ho  and  Probstein,  is  qualitatively  confirmed  in  that  some 

measured  values  of  (A/Rb)/(^/Rb)inviscid  less  than  unity.  The 
constant- density  solution  of  Probstein  and  Kemp  should  approximate 
more  closely  to  the  uncooled  model  condition  than  to  the  cold  wall  case; 

Fig.  6  shows  this  to  be  true.  It  is  important  to  remember  that  the  re¬ 
sults  obtained  by  Ho  and  Probstein  and  Probstein  and  Kernp  apply  to  air, 
whereas  the  present  data  and  the  results  of  Levinsky  and  Yoshihara 
apply  to  argon. 

It  is  considered  pertinent  here  to  emphasize  that  the  value  A  quoted 
herein  is  the  distance  from  the  face  of  the  body  to  the  leading  edge  of  the 
shock.  For  high  Reynolds  numbers  the  flow  in  front  of  a  body  is  divided 
into  three  regions;  the  shock  layer,  the  boundary  layei',  and  the  distance 
from  the  body  to  the  rear  face  of  the  shock  (A).  When  the  boundary  and 
shock  layers  merge,  these  three  divisions  no  longer  have  any  significance 
as  separate  entities.  For  this  reason,  the  shock  detachment  distance.  A, 
is  taken  to  be  the  distance  from  the  face  of  the  body  to  the  shock  leading 
edge. 

4.1,3  Flat-Faced  Body  Results 

Some  limited  tests  wore  made  with  two  water-cooled,  flat-faced 
probes  at  the  same  flow  conditions  as  for  the  water-cooled  spheres.  The 
results  of  the  tests  are  shov/n  in  Fig.  3.  No  attempt  has  been  made  to 
draw  a  single  curve  through  these  results  since  the  sphere  data  imply 
that  there  will  be  a  variation  with  Mach  number,  which  the  limited  nature 
of  these  tests  of  flat-faced  bodies  cannot  fully  specify.  , 

The  results  of  this  test  are  interesting  when  they  are  considered  in 
conjunction  with  the  sphere  data.  First  of  all,  the  shock  detachment 
distance  is  more  than  twice  the  value  obtained  for  a  sphere  at  the  same 
Mach  number  and  Reynolds  nuiuber.®  This  implies  that,  provided  the 
boundary  layer  increases  in  the  same  manner  as  for  a  sphere,  the  boundary 
and  shock  layers  will  not  merge  until  a  lower  Reynolds  number  than  that 
at  which  merging  occur, s  for  a  sphere  is  reached.  Thus,  incipient  and 
merged-layer  phenomena  would  not  be  expected  to  occur  with  flat-face 
bodies  until  a  lower  Reynolds  number  than  that  at  which  they  occur  for 
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spheres  is  reached.  In  connection  with  this,  Sherman's  (Ref.  16)  work 
Oii  impact-pressuie  probes  shows  that  the  rapid  increase  in  impact  pres¬ 
sure  of  Reynolds  number  occurs  at  higher  Reynolds  number  for  source¬ 
shaped  probes  than  for  flat-faced  probes.  This  is  consistent  with  the 
idea  that  merged-layer  characteristics  occur  at  higher  Reynolds  numiber 
for  source-shaped  than  for  flat-faced  probes.  ® 

It  has  been  shown  that  shock  detachment  distance  is  a  function  of 
wall  temperature,  among  other  things.  In  Ref.  17,  a  study  was  made 
of  flat-faced,  impact-pressure  probes  where  the  wall  to  stagnation 
temperature  ratio  w'as  approximately  0.  3.  It  was  found  that  these  results 
did  not  compare  very  well  with  some  data  obtained  for  a  wall  to  stagnation 
temperature  ratio  of  unity.  The  present  investigation  would  seem  to  indi¬ 
cate  that  this  difference  could  be  at  least  partly  accounted  for  by  the  dif¬ 
ference  in  shock  detachment  distance  for  the  two  cases.  Van  Dj'ke  (Ref,  1) 
has  shown  that  shock  dftachrnent  distance  is  primarily  a  function  of  the 
density  ratio  across  a  normal  shock  wave.  For  a  sphere,  the  effect  of 
different  specific  heat  ratios  is  small.  For  the  density  ratios  of  the  pres¬ 
ent  tests,  the  error  involve<;^in  the  shock  detachment  distance  by  taking 
;■'  =  1.4  inslead  of  i.  667  is  on  the  order  of  3  percent.  Assuming  that  the 
same  is  true  for  flat-faced  bodies,  then  a  good  approximation  to  the 
inviscid  value  of  A/R^  for  the  flat-faced  bodies  can  be  obtained  from  ex¬ 
perimental  data  for  flat-faced  bodies,  and  this  is  shown  in  Fig,  8.  The 
inviscid  values  of  A/Rb  for  the  density  ratios  of  the  present  tests  are 
shown  in  Fig.  3.  As  for  the  spheres,  there  is  an  indication  that  for  the 
Reynolds  numbers  of  the  present  tests  the  measured  value  of  A/Rb  is  less 
than  the  inviscid  value  through  an  intermediate  range  of  Reynolds  numbers. 

4.2  SHOCK  SHAPE 
4.2.1  Sphere  Results 

The  results  of  the  analysis  to  determine  the  shock  shape  are  shown 
in  Figs.  9  and  10.  For  the  range  of  Macii  numbei-  and  Reynolds  number 
of  these  tests,  the  nondimensionalized  shock  shape  is  a  function  of  Mach 
number  only,  after  the  effect  of  Reynolds  number  on  A  is  accounted  for. 

In  order  to  account  for  the  Reynolds  number  effect  on  the  shock  detach¬ 
ment  distance.  A,  the  origin  of  the  ordinates  is  taken  to  be  at  the  apex 
of  the  shock-wave  profile. 

In  Ref.  18,  an  analysis  of  experimentally  determined  shock  shapes 
has  shown  that  these  shapes  could  be  represented  by  an  equation  of  hyper¬ 
bolic  form: 


O 
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withthe  apex  at  the  origin  and  with  A  being  the  slope  of  the  asymptotic 
shock  angle  far  from  the  body.  When  n  =  1/2,  this  equation  is  exactly 
hyperbolic.  From  the  above-mentioned  shock-shape  analysis,  values 
were  assigned  to  A,  n,  and  a  and  were  shown  to  be  functions  of  free- 
stream  Mach  number  and  the  nose  wave  drag  coefficient,  calculated  by 
using  the  modified  Newtonian  approximation  and  the  specific  heat  ratio 
of  the  gas.  On  comparing  the  shock  shapes  predicted  by  this  method 
with  those  calculated  by  an  empirical  relationship  by  Love  (Ref.  19), 
it  was  found  that  Love's  method  gave  a  better  prediction  for  Mach  num¬ 
bers  less  than  four.  The  empirical  relationships  derived  in  Refs.  18 
and  19  are  for  air  only,  and  neither  of  them  can  be  simply  modified  to 
account  for  a  gas  with  a  specific  heat  ratio  other  than  1.4.  It  has  been 
shown  in  Ref.  6  that  the  blast  analogy  is  a  useful  means  for  predicting 
some  of  the  features  of  hypersonic  flow  over  blunt  bodies.  In  fact,  it 
is  possible  with  this  method  to  predict  the  relative  effects  of,  say,  a 
change  in  specific  heat  ratio  very  well,  while  being  unable  to  predict 
a  precise  value  of  the  quantity  under  investigation.  In  terms  of  the 
blast  analogy,  the  shock  shape  in  front  of  a  body  can  be  written  as: 


This  is  a  useful  form  when  a  comparison  of  the  effects  of  different 
specific  heat  ratios  on  shock  shape  is  to  be  made.  When  the  relevant 
values  of  the  specific  heat  ratio  for  air  and  argon  are  substituted  in  this 
equation  at  constant  value  of  x^/D,  the  ratio  (y/D)aj.gon/(y/^)air  ""  ^6. 

Some  experimental  results  for  helium  (Refs.  20  and  21),  which  has  the 
same  specific  heat  ratio  as  argon  over  a  Mach  number  range  from  15  to 
40,  show  that  this  ratio  is  1.17. 

For  the  present  Mach  number  range  of  4.  22  to  9.  21,  the  method 
given  in  Ref.  19  has  been  used  to  determine  the  shock  shape  for  spherical 
bodies  in  air.  The  data  for  Mach  numbers  of  8.  61  and  9.  21  were  ob¬ 
tained  by  extrapolating  the  results  contained  in  Ref.  19.  It  was  considered 
that  this  was  preferable  to  using  another  method,  say  Ref.  18,  to  calculate 
the  shock  shape  at  these  high  Mach  numbers.  The  resulting  shock  shapes 
are  shown  plotted  in  Figs.  9  and  10.  The  ratio  (y/ D)g^gQ^/(y derived 
from  the  experimental  argon  data  and  the  calculated  air  data  is  shown 
plotted  against  Mach  number  in  Fig.  11.  At  the  high  Mach  number  end  of 
this  curve,  this  ratio  agrees  very  well  with  the  blast  analogy.  Intuitively, 
it  would  seem  that  the  ratio  must  decrease  with  decreasing  Mach  number. 
No  reason  can  be  found  for  the  abrupt  change  in  slope  of  the  curve  at  a 
O  Mach  number  of  approximately  six. 
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4.2.2  Flat-Faced  Results  ^ 

The  measured  shock  shapes  for  the  flat-faced  bodies  are  shown  in 
Fig.  12.  The  shock  shapes  shown  for  air  were  derived  from  some  ex¬ 
perimental  data  for  flat-faced  bodies  contained  in  Refs.  19  and  22. 
Equation  (7)  indicates  that  the  ratio  (y/ D)gj.gQn/(y/D)air  should  be  the 
same  for  spheres  and  for  flat-faced  bodies.  The  dashed  lines  in  Fig.  12 
were  obtained  by  using  the  experimentally  determined  curves  for  air  and 
multiplying  them  by  the  experimentally  determined  ratio  in  Fig.  11.  It 
will  be  seen  that  the  dashed  curves  fit  the  experimentally  measured 
shock  shapes  for  flat-faced  bodies,  which  provides  additional  confirma¬ 
tion  of  the  form  of  the  curve  shown  in  Fig.  1 1 . 

4.2.3  Correlotion  of  Shock  Shapes 


In  terms  of  the  blast  analogy  the  parameter  governing  the  form  of 
the  shock-wave  shape  at  high  Mach  numbers  is  given  by  Eq.  (7)  as 

for  a  constant  value  of  y.  For  Mach  numbers  greater  than 
4.  0  the  inviscid  drag  coefficient  of  both  flat-faced  bodies  and  spheres 
is  approximately  constant.  Therefore,  to  correlate  the  experimental 
shock  shape  data  for  spheres  and  flat-faced  bodies,  it  should  be  suffi¬ 
cient  to  plot  x/D  versus  y/D  on  log-log  paper.  Figure  13  shows  the 
result  of  plotting  data  obtained  from  the  faired  curves  of  Figs.  9,  10, 
and  12,  for  both  argon  and  air.  The  fact  that  a  series  of  parallel  straight 
lines  can  be  drawn  through  the  sphere  data  indicates  that  a  constant  value 
can  be  given  to  the  exponent  of  (xij/D)  in  Eq.  (7).  However,  (I'coCd/^Jq) 
will  be  a  function  of  Mach  number  for  ^  6,  a  fact  which  is  not  pre¬ 
dicted  by  the  blast  analogy  which  applies  to  hypersonic  flows.  A  straight 
line  cannot  be  drawn  through  the  flat -faced  body  data  for  values  of 
x/D  <  1.  0,  which  means  that  a  constant  value  cannot  be  assigned  to  the 
exponent  in  Eq.  (7).  There  are  insufficient  data  to  say  definitely  whether 
there  is  a  linear  relationship  between  x/D  and  y/D  for  values  of  y/D  >1.  0, 
though  there  are  indications  to  suggest  that  this  is  so. 

4.3  SHOCK-WAVE  THICKNESi 

• 

In  Section  4.  1. 2,  a  theory  derived  by  Ziering,  et  al.  (Ref.  11),  for 
monatomic  gases  was  used  to  estimate  the  thickness  of  the  shock  wave 
for  the  present  test.s.  This  estimate  was  chosen  because  it  appeared  to 
be  in  reasonable  agreement  with  the  experimental  data  for  Mach  numbers 
up  to  5.  A  comparison  of  some  of  the  theoretical  estimates  of  shock-wave 
thickness  for  monatomic  gases  is  shown  in  Fig.  14  together  with  the  ex- 
perim.ental  data  published  by  other  investigators. 

In  the  calci  'ation  of  the  shock-wave  thickness,  two  approaches  have 
been  used.  In  one,  the  Navier-Stokes  equations  have  been  used  to  calculate 
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shock-wave  profiles,  and  curve  6  in  Fig.  14  shows  a  solution  by 
Schwartz,  as  reported  by  Hornig  in  Ref.  23.  The  other  solutions  are 
developments  of  a  method  suggested  by  Mott-Smith  (Ref.  24)  for  the 
solution  of  the  Boltzmann  equation.  Curves  1-5  in  Fig.  14  (Refs.  11, 

23,  25-27)  give  the  results  of  calculations  of  this  form.  For  Mach 
numbers  less  than  approximately  two,  all  of  the  calculations  are  in 
reasonable  agreement  with  the  available  experimental  data.  However, 
for  higher  Mach  numbers,  the  degree  of  agreement  is  poor.  At  a  Mach 
number  of  six,  the  solution  proposed  by  Koga  gives  a  value  of  the  shock 
thickness  five  times  larger  than  that  predicted  by  the  Navier-Stokes 
equations,  with  the  other  solutions  lying  within  this  range. 

O 

By  far  the  most  consistent  and  extensive  set  of  experimental  data 
to  date  is  that  by  Linzer,  reported  by  Hornig  in  Ref.  23.  These  values 
of  shock-wave  thickness  have  been  derived  from  measurements  of  the 
optical  reflectivity  of  the  shock  wave.  Although  most  of  the  other  test 
data  reported  in  Ref.  28  are  in  good  agreement  with  the  above  data  at 
low  Mach  numbers,  there  is  some  scatter  in  the  data  at  the  higher 
Mach  numbers.  In  spite  of  this  scatter,  it  seems  fair  to  conclude  from 
the  experimental  data  that  the  shock-wave  thickness  has  a  minimum 
value  at  M  -  3.  5  and  thereafter  increases  with  increase  in  Mach  number. 
This  form  of  variation  is  the  same  as  that  predicted  by  Ziering,  et  al. 

(Ref.  11),  although  the  actual  values  of  the  shock  thickness  are  not  the 
same. 

• 

jjt  is  important  to  remember  that  the  shock- wave  thickness  con¬ 
sidered  here  is  the  thickness  of  the  density  profile.  Furthermore,  the 
thicknesses  quoted  are  the  maximum  slope  thicknesses  of  this  profile, 
as  illustrated  in  Fig.  15.  Sherman  (Ref.  29)  has  shown  that  there  is 
considerable  difference  between  the  total  shock-wave  thickness  and  the 
maximum  slope  thickness  in  the  low-density  flow  regime.  In  fact,  when 
consideration  is  being  given  to  the  merging  of  shock  and  boundary  layers 
in  front  of  a  blunt  body,  it  would  seem  that  the  total  thickness  is  more 
significant  than  the  maximum  slope  thickness. 

• 

Consider  the  structure  of  the  shock  wave  at  locations  removed  from 
the  spheres  and  flat-faced  bodies  of  the  present  tests.  If  it  is  assumed 
that  the  film  density  of  the  photographic  negatives  is  proportional  to  the 
density  in  the  gas  flow,  then  the  measurement  of  the  film  density  witli^he 
photo- densitometer  should  give  the  form  of  the  density  variation  through¬ 
out  the  flow  regime  under  investigation.  As  stated  earlier  for  all  the 
flow  conditions  investigated  herein,  there  is  an  axial  Mach  number  gra¬ 
dient,  which  in  turn  implies  a  density  gradient.  Measurement  of  the  film 
density  upstream  of  the  body  shock  indicates  that  the  film  density  decreases 
as  the  shock  is  approached,  i.  e. .  in  passing  downstream.  Furthermore, 
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the  rate  of  decrease  of  film  density  decrea.ses  as  the  region  of  flow 
under  investigation  moves  farther  and  farther  downstream  from  the 
nozzle  exit.  This  implies  a  decrease  of  axial  Mach  number  gradient 
as  the  distance  from  the  nozzle  exit  plane  increases.  This  is  consist¬ 
ent  with  the  results  of  the  nozzle  calibration  studies  which  show  just 
such  a  decreasing  gradient. 

The  film  density  of  a  large  number  of  the  negatives  was  recorded 
on  the  strip  chart  recorder  attached  to  the  photo-densitometer  for  a 
series  of  stations  along  the  bow  shock,  but  sufficiently  far  from  the 
bodies  to  avoid  boundary-layer  influences.  As  the  shock  was  approached 
from  upstream,  the  density  decreased,  then  underwent  a  gradual  increase 
at  the  front  of  the  shock,  then  a  rapid  rise,  then  a  more  gradual  rise, 
then  became  almost  constant,  and  finally  the  density  started  to  decrease. 

An  example  is  shown  in  Fig.  16.  Consider  the  density  variation  along 
section  AA  through  the  three-dimensional  shock  shown  in  Fig.  17.  From 
a-b,  the  density  decreases  because  of  the  axial  Mach  number  gradient; 
from  b-c-d,  an  increase  occurs  in  density  through  the  shock  to  a  maxi¬ 
mum  value  at  d  at  the  rear  face  of  the  shock;  and  from  d-e,  a  decrease  in 
density  occurs  because  of  the  flow  expansion  behind  the  shock.  The  form 
of  this  density  variation  is  the  same  as  that  recorded  by  the  photo¬ 
densitometer  in  Fig.  16.  Now,  if  the  radiation  of  light  from  the  shock 
wave  is  considered,  it  will  be  seen  that  the  maximum  section  length 
through  the  shock  wave  occurs  at  section  d.  This  indicates  that  the  points 
of  maximum  film  density  and  flow  density  can  occur  at  the  same  points. 

By  adjusting  the  rate  at  which  the  carriage  carrying  the  negative  is  passed 
by  the  viewing  slit  of  the  densitometer,  the  shock  thickness  recorded  on 
the  strip  chart  recorder  can  extend  over  several  inches.  Both  the  maxi¬ 
mum  slope  and  total  thickness  of  the  shock  waves  were  measured  from 
this  recorded  data.  These  thicknesses  were  measured  normal  to  the 
tangent  to  the  shock  wave.  Then,  taking  the  free- stream  mean  free  path 
and  Mach  number  resolved  in  a  direction  normal  to  the  shock  wave,  the 
ratio  of  the  maximum  slope  shock-wave  thickness  to  the  free- stream 
mean  free  path  was  calculated  and  is  shown  plotted  in  Fig.  18,  The  scatter 
in  the  dhta  is  on  the  order  of  ±10  percent.  This  shows  that  as  the  Mach 
number  increases  above  roughly  three,  the  shock-wave  thickness  increases, 
which  is  in  agreement  with  the  form  predicted  in  Ref.  11,  An  attempt  has 
been  made  to  estimate  the  total  thickness  even  though  it  was  difficult  to 
define  exactly  where  the  edges  of  the  shock  wave  were.  Figure  18  shows 
that  the  total  thickness  is  approximately  50  percent  greater  than  the  maxi¬ 
mum  slope  thickness.  Sherman  {Ref.  29)  has  made  a  study  of  the  shock- 
wave  thickness  using  a  fine  wire  probe  (which  was  essentially  a  frec- 
molecule  flow  probe)  in  helium  and  air.  The  results  of  this  study  in  the 
Mach  number  range  from  1.  72  to  3,  91  show  that  the  total  thickness  is 
always  approximately  double  the  maximum  sl^e  thickness.  Also  of  interest 
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is  the  fact  that  in  this  study  the  wire  temperature  profile  is  not  sym¬ 
metrical  about  the  mean  temperature  value.  Some  of  the  film  density 
profiles  obtained  in  the  present  tests  were  also  of  this  form.  Hornig 
(Ref.  23)  indicates  that  theory  predicts  that  this  asymmetry  is  a  prop¬ 
erty  of  strong  shocks.  It  would  seem  from  the  present  results  and 
those  of  Sherman  that  this  property  of  strong  shock  waves  is  confirmed. 
As  mentioned  by  Hornig,  such  a  property  of  shock  waves  poses  a  prob¬ 
lem  in  the  optical  reflectivity  method  of  shock-wave  thickness  measure¬ 
ments,  since  the  basis  for  the  method  is  the  assumption  that  the  density 
profile  through  the  shock  is  symmetrical. 

9 


4.4  SPHERE  PRESSURE  DISTRIBUTIONS 


Surface  pressure  distributions  in  the  Mach  number  rangg  from  4.  4 
to  13.  5  are  shown  in  Fig.  19  and  Table  i.  A  compai'ison  is  made  with 
the  modified  Newtonian  approximation; 


This  equation  can  be  rewritten  in  terms  of  surface  pressure  and  the 
pressure  at  the  stagnation  point: 


-  cos^  0  4  -Hsi—  sin^  0  (g) 

Po 

For  the  complete  Mach  and  Reynolds  number  range  of  these  tests, 
the  modified  Newtonian  approximation  gives  a  good  prediction  of  the 
surface  pressure  distribution  up  to  an  angle  of  60  deg.  Above  60  deg, 
there  is  a  tendency  for  the  modified  Newtonian  approximation  to  under¬ 
estimate  the  surface  pressure  distribution.  However’,  it  must  be 
emphasized  that  the  degree  of  disagreement  between  theory  and  experi¬ 
ment  is  small  even  at  surface  angles  greater  than  60  deg. 

Talbot,  et  al.  (Ref.  30),  have  shown  for  hemispherically  blunted 
cones  in  an  unheated,  lov/-donsity  flow  that  the  modified  Newtonian 
approximation  gives  good  results  over  the  hemispherical  part  of  the 
body.  Hayes  and  Probstein  (Ref.  31),  in  a  comparison  of  the  merits  of 
the  various  methods  of  predicting  the  surface  pressure  distribution  over 
spheres,  emphasize  that  the  m.odified  Newtonian  approximation  is  strictly 
an  empirical  formula,  A  more  rational  estimate  is  the  Newton-Busemann 
law  which  takes  into  account  centrifugal  effects: 


This  law,  though  more  rational  in  nature  than  the  modified  Newtonian 
approximation,  gives  poor  agreement  with  experiment.  Hayes  and 
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Probstein  (Ref.  31)  consider  that  the  good  results  obtained  with  the 
modified  Newtonian  approximation  are  fortuitous  and  that  the  centrif¬ 
ugal  pressure  difference  across  the  shock  layer  is  in  many  instances 
offset  by  the  effect  of  the  difference  in  the  shock  and  body  angles. 

Perhaps  one  of  the  more  interesting  results  from  this  experiment 
is  that  the  modified  Newtonian  theory  is  valid  in  a  region  where  the 
shock  and  boundary  layers  merge.  The  only  other  known  data  in  the 
low-density  regime  where  the  shock  and  boundary  layers  had  not 
merged  are  contained  in  Ref.  30. 


It  is  of  interest  to  consider  the  form  of  the  surface  pressure  dis¬ 
tribution  in  the  free-molecule  flow  regime  in  order  to  determine  how 
the  form  of  the  distribution  changes  in  going  from  the  continuum  to 
free-molecule  regimes.  The  surface  pressure  distribution  of  a  sphere 
in  the  free-molecule  flow  regime,  given  in  Ref.  32,  is 


C  09 


’<9 


+ 


V  s 


cos  6  (1  +  erf  S  cos  6) 


I 


(11) 


If  the  pressure  at  the  stagnation  point  is  written  as  Po',  then  the  pres 
sure  at  any  point  on  the  sphere  can  be  written  as 


P*  ^  [g-S^  cos^e  ^  s  cos  0  (  1  +  erf  S  cos  9) 
P„  '  [  c“®*  4  /jT  .S  (  1  4  erf  S)  j 


for  S  cos  0  >  2.0 

then  erf  .S  cos  0  -<•  1.0 

then  Eq.  (16)  becomes 


(p7) 


when  .S  cos  e  >  2 


(12) 


(13) 


For  the  present  series  of  tests,  this  would  imply  that  if  any  free- 
molecule  flow  effects  were  present,  the  measured  pressures  would  be 
larger  than  those  predicted  by  the  modified  Newtonian  approximation. 
Since  the  measured  pressures  are  not  larger  than  those  predicted  by 
the  modified  Newtonian  approximation,  then  the  spheres  are  not  being 
affected  by  any  free-molecule  flow  effects.  This  is  not  a  surprising 
result  when  it  is  considered  that  the  largest  Knudsen  number,  based  on 
sphere  radius,  of  the  present  tests  is  approximately  0,  07  for  a  billiard- 
ball  molecular  model. 
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5.0  CONCLUSIONS 


Based  on  an  experimental  investigation  to  determine  the  pressure 
distribution,  shock  shape,  and  shock  detachment  distance  for  spheres 
and  the  latter  two  characteristics  for  flat-faced  bodies  in  heated  argon 
flow,  the  following  conclusions  are  made; 

The  shock  detachment  distance  in  front  of  a  body  in  low- density, 
hypervelocity  flow  is  a  function  of  Mach  number,  wall-to-stagnation 
temperature  ratio,  and  Reynolds  number  and  body  shape.  Evidence 
from  the  results  of  the  present  series  of  experiments  suggests  that  the 
shock  wave  and  boundary  layers  in  front  of  the  bodies  were  merged. 

For  complete  range  of  Ma^h  and  Reynolds  numbers  encounterec^ 
in  this  experiment,  the  nondimensionalized  shock  shape  was  a  function 
of  Mach  number  and  body  shape  only.  Furthermore,  the  difference  be¬ 
tween  the  shock  shapes  in  air  and  argon  at  high  Mach  numbers  was  in 
very  good  agreement  with  that  predicted  by  the  blast  analogy. 

s 

Even  in  the  flow  regime  where  the  shock  and  boundary  layers 
were  merged,  the  modified  Newtonian  approximation  gave  good  results 
for  the  surface  pressure  distribution  over  spheres  for  6  ■'  60  deg. 

Measurements  of  shock-wave  thickness  in  the  low-density  regime 
were  in  reasonable  agreement  with  results  at  higher  densities  when 
ratioed  to  mean  free  path.  Also,  the  form  of  the  variation  of  shock- 
wave  thickness  with  Mach  number  agreed  with  that  predicted  by  Ziering, 
et  al,  (Ref.  11). 
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Argon  -  Schwartz  and  Eckerman  (Ref.  14) 
Helium  -  Vas,  et  al .  (Ref.  15) 


Fig.  6  Shock  Detachment  Distance  for  a  Sphere  in  a  Monatomic  Gas  at  High  Reynolds  Number 
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Fig.  7  Comparison  of  Theoretical  and  Experimental  Shock  Detachment  Distances 
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Fig.  10  Continued 
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o.  M  =  4.22,  To  =  2780°K 
Fig.  12  Shock  Shape  for  Water-Cooled,  Flot-Faced  Body 
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Fig.  14  Vorlatlon  of  ^hock-Wave  THtckness  with  Moch  Number  for  Monatomic  Goses 
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Fig.  15  Shock-Wave  Density  Profile 
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Arnold  Engineering  Development  Center  1.  Shock  waves 
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